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ABSTRACT 



We model the nucleosynthesis during the thermal pulse phase of a rotat- 
ing, solar metallicity Asymptotic Giant Branch (AGB) star of 3M , which was 
evolved from a main sequence model rotating with 250km/s at the stellar equa- 
tor. Rotationally induced mixing during the thermal pulses produces a layer 
(~ 2 • 1O~ 5 M ) on top of the CO-core where large amounts of protons and 12 C 
co-exist. With a post-processing nucleosynthesis and mixing code, we follow the 
abundance evolution in this layer, in particular that of the neutron source 13 C 
and of the neutron poison 14 N. In our AGB model mixing persists during the 
entire interpulse phase due to the steep angular velocity gradient at the core- 
envelope interface, thereby spreading 14 N over the entire 13 C-rich part of the 
layer. We follow the neutron production during the interpulse phase, and find 
a resulting maximum neutron exposure of r max = 0.04mbarn _1 , which is too 
small to produce any significant s-process. In parametric models, we then inves- 
tigate the combined effects of diffusive overshooting from the convective envelope 
and rotationally induced mixing. Just adding the overshooting and leaving the 
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rotational mixing unchanged results also in a small maximum neutron expo- 
sure (O.OSmbarn" 1 ). Models with overshoot and weaker interpulse mixing — as 
perhaps expected from more slowly rotating stars — yield larger neutron expo- 
sures. A model with overshooting without any interpulse mixing obtained up to 
0.72mbarn _1 , which is larger than required by observations. We conclude that 
the incorporation of rotationally induce mixing processes has important conse- 
quences for the production of heavy elements in AGB stars. While through a 
distribution of initial rotation rates, it may lead to a natural spread in the neu- 
tron exposures obtained in AGB stars of a given mass in general — as appears 
to be required by observations — it may moderate the large neutron exposures 
found in models with diffusive overshoot in particular. Our results suggest that 
both processes, diffusive overshoot and rotational mixing, may be required to 
obtain a consistent description of the s-process in AGB stars which fulfils all ob- 
servational constraints. Finally, we find that mixing due to rotation within our 
current framework does increase the production of 15 N in the partial mixing zone. 
However, this increase is not large enough to boost the production of fluorine to 
the level required by observations. 

Subject headings: stars: AGB and post-AGB — stars: evolution — nuclear 
reactions, nucleosynthesis, abundances — stars: rotation — stars: interiors 

1. Introduction 

Trans-iron elements are mainly made by neutron-capture reactions on 56 Fe seed nuclei. 
Two processes have been been distinguished according to the neutron density at the pro- 
duction site. In the case of the r-process the n-densities are high (N n > 10 20 cm~ 3 ), and the 
time scale of successive n-capture reactions on heavy isotopes is faster than the /3-decay time 
scale. Such a sudden high-density neutron burst creates isotopes far away from the valley 
of /^-stability in the chart of nuclides, which successively decay back to the stable isotopes. 
In contrast the s-process is characterized by lower neutron densities (N n ;$ 10 10 cm~ 3 ). Neu- 
tron captures are generally followed by /3-decays since unstable isotopes on the s-process 
path have typical life times of the order of hours. In some cases, however, the unstable 
isotopes involved have longer life times and, depending on the neutron density and temper- 
ature conditions, branchings can be open on the s-process path leading to the production 
of neutron-rich isotopes (see Clayton 1968, for an introduction to the s-process). 

In Asymptotic Giant Branch (AGB) stars recurrent He-shell flashes (thermal pulses, TP) 
drive a convective zone that temporarily covers the whole region between the H-burning and 
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the He-burning shells (intershell). Here, partial He burning produces a high mass fraction 
(> 0.25) of 12 C and the chain 14 N(a,7) 18 F(/3+z/) 18 (a,7) 22 Ne starting on the abundant 14 N 
from the H-burning ashes produces a relatively large amount of 22 Ne (mass fraction ~ 0.02). 
The 22 Ne(a, n) 25 Mg reaction as neutron source for the s-process was suggested by Cameron 
(1960). Temperatures above 3-10 8 K are required for that reaction to release a significant 
amount of neutrons. Stellar models showed that such high temperatures are achieved in 
intermediate-mass (M Z ams/M > 4) AGB stars, which were hence proposed as the main 
site for the production of s-process elements belonging to the solar main component, i.e. 90 
< A <204 (Iben 1975; Truran & Iben 1977). However, the neutron density produced by 22 Ne 
burning in thermal pulses is rather high (above 10 u cm~ 3 for T = 3.5 • 10 8 K). This leads 
to excesses in the neutron-rich nuclides produced by branchings, for example 86 Kr, 87 Rb, 
and 96 Zr (see e.g. Despain 1980), in contrast with the great majority of the observations 
of s-pro cess-enhanced stars such as MS, S and C stars. In S and C stars the Rb/Sr ratio 
is typically much lower than solar (Lambert et al. 1995; Abia et al. 2001), indicating a 
low neutron density at the s-process site. Also the Rb abundances observed in 10 AGB 
members of the massive Galactic globular cluster uj Centauri indicate a low neutron density 
for the s-process (Smith et al. 2000). Lambert et al. (1995) reported the zirconium isotopic 
abundance obtained by spectroscopic observations of the ZrO bandheads in M, MS and S 
stars and found no evidence of an excess of the neutron-rich 96 Zr which can be produced in 
great amount by the s-process when the neutron density exceeds ~ 5 • 10 8 n cm -3 . Another 
problem is that the activation of the 22 Ne(a, n) 25 Mg is expected to produce an excess of 
25 Mg in stars enriched in s-process elements. Instead these stars typically have magnesium 
isotopic abundances in solar proportion (see e.g. Smith & Lambert 1986; McWilliam & 
Lambert 1988). Also other types of observations tend to exclude AGB stars of intermediate 
mass as the main s-process site. Observations mainly show that MS, S and C stars have 
low luminosity (Frogel et al. 1990) and hence low mass. Feast (1989) performed a study of 
the kinematics of peculiar red giants including S, SC, and C stars. On the basis of 183 S - 
SC stars, and 463 C stars he estimated their mean mass to be 1.3 and 1.6 M respectively, 
although this estimate needs to be improved. In summary, the observational evidence and 
the current state of AGB evolution models suggest that the major nuclear production site 
of the s-process are low-mass AGB stars. 

In low-mass AGB stars the temperature in the intershell is not high enough to burn a 
significant amount of 22 Ne. The 13 C(o;,n) 16 neutron source reaction, which was suggested 
by Greenstein (1954) and Cameron (1955) and is activated at lower temperatures (~ 0.9 • 10 8 
K), is expected to play the major role (Iben & Renzini 1982; Gallino et al. 1988; Kappeler 
et al. 1990). However, an amount of 13 C higher than that present in the H-burning ashes 
is needed to reproduce the observed enhancements of heavy elements. In order to form 
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a sufficient amount of 13 C it is hence speculated that some protons mix into the 12 C-rich 
intershell (see Busso et al. 1999, for a general review of the s-process in AGB stars). 

In recent years a picture of the s-process based on these results has emerged (Gallino 
et al. 1998) and is summarized in Figure 1. The He-flash convection zone homogenizes 
the intershell region and 12 C produced in the He-burning shell is mixed up to just below 
the location of the H-burning shell. The dashed line in Figure 1 indicates that, after the 
convective pulse is extinguished, the convectively unstable envelope may extend down into 
H-free layers of the intershell region. This phenomenon allows processed intershell material 
to be carried into the envelope and hence to the stellar surface (third dredge-up). At the 
end of third dredge-up a layer is created where the H-rich envelope directly neighbours the 
12 C-rich intershell. This layer is a favourable region for the formation of a zone where 12 C 
and protons are partially mixed. As the temperature increases in the region a pocket of 13 C 
forms by proton captures on 12 C. Subsequently, the 13 C serves as a neutron donor via the 
reaction 13 C (a, n) 16 which is activated during the following interpulse phase at T ~ 9- 10 7 K 
so that neutrons are released under radiative conditions (Straniero et al. 1995). Typically the 
s-process occurs on a time scale of several tens of thousands years and before the onset of the 
following thermal pulse. In the convective pulse the 22 Ne neutron source is only marginally 
activated. 

It is most reasonable to assume that the H/ 12 C-ratio in the partial mixing zone varies 
somehow continuously from a few hundreds in the envelope to zero in the intershell. Then, 
one finds in the top layers of the partial mixing zone another pocket made of 14 N, which 
forms where the H/ 12 C-ratio is larger than where the 13 C pocket forms. Without further 
mixing the partial mixing layer is strictly stratified during the interpulse period: as shown in 
Figure 1 the two pockets coexist without much interaction. During the interpulse phase the 
temperature does not reach values required for the 14 N(a,7) 18 F reaction. The 14 N pocket 
is engulfed into the following He-flash convection zone where it might slightly contribute to 
the production of 22 Ne. 

Stellar models that use a standard treatment of convective mixing can not produce the 
13 C pocket because extra mixing processes are required to allow a small amount of protons 
to enter the 12 C-rich region. Gallino et al. (1998, and following studies by that group) and 
Goriely & Mowlavi (2000) have studied the s-process by assuming a certain proton profile 
extending into the 12 C-rich region without relating explicitly to a specific physical process. 
In these studies it was implicitly assumed that any subsequent alterations of the abundances 
in the partial mixing zone are due to nucleosynthesis only and no mixing takes place during 
the interpulse phase after the initial formation of the partial mixing zone. In this way it was 
possible to develop a fairly good understanding of the s-process nucleosynthesis in the 13 C 
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pocket. For example, it was shown that many observations are reproduced with models in 
which the neutron exposure in the 13 C pocket is up to ~ 0.4mbarn _1 at solar metallicity. In 
order to account for the observed s-process overabundances the partial mixing zone needs to 
have a mass of 10~ 4 - 10~ 3 M Q . These s-process models with an assumed H-profile for the 
partial mixing zone can account for many of the overall observed properties of the s-process. 

However, to explain the observed s-process signature in AGB stars as a function of 
metallicity, the s-process model described above requires that stars of the same mass and 
metallicity have different neutron exposures due to the 13 C neutron source during the AGB 
interpulse phase (Busso et al. 2001). A mixing process for the formation and evolution 
of this neutron source that allows for some spread in the efficiency of producing neutrons 
seems to be necessary. Van Winckel & Reyniers (2000) arrive at the same conclusion from 
observations of weakly metal-poor post-AGB stars with 21/xm feature. In these objects the 
s-process element signatures are easier to measure than in the cool progenitor AGB stars. 
Also the spread in the Pb abundance observed in very low-metallicity stars (for example 
Van Eck et al. 2003) could fit into the current s-process model only if there is a spread 
in the number of neutrons available in the pocket for the s-process. Such a spread in the 
efficiency of the 13 C neutron source is also required to explain the measurements of isotopic 
ratios in single pre-solar silicon carbide grains (SiC) recovered from pristine carbonaceous 
meteorites (Zinner 1998). The majority of these grains are believed to have formed in the 
circumstellar dusty envelopes surrounding carbon stars. The main evidence for this comes 
from measurements of aggregates of SiC that have shown a very strong signature of the 
s-process in the heavy elements present in trace (see e.g. Gallino et al. 1997). Nicolussi et al. 
(1997, 1998a,b) and Savina et al. (2003) performed measurements of the composition of Sr, 
Zr, Mo and Ba in single SiC grains. Isotopic ratios that are sensitive to the efficiency of 
the 13 C neutron source have been found to show a large spread within the single grain data. 
This can be explained if SiC grains were formed in a multiplicity of low mass AGB stars with 
a range of neutron exposures in the 13 C pocket. These isotopic ratios are those involving 
nuclei with very low neutron capture cross section such as nuclei with number of neutrons 
equal or near to nuclear magic numbers, i.e. 88 Sr/ 86 Sr, 90 > 91 > 92 Zr/ 94 Zr and 138 Ba/ 136 Ba. A 
detailed analysis of SiC grain data indicated that a spread of the order 5 is necessary in the 
neutron exposure for a given mass and metallicity (Lugaro et al. 2003a). These new insights 
are specifically important for constraining the physical processes which are responsible for 
the partial mixing between envelope and core, and thereby lead to the formation of the 13 C 
neutron source. Moreover, neutron capture elements are in general becoming an increasingly 
important target of stellar observations, for example of IR-observations in planetary nebulae 
(Dinerstein 2001). 

The open problem addressed in this paper is the role of rotational mixing. Rotation is an 
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effect that has to be taken into account when studying stellar evolution and nucleosynthesis. 
Most F stars, which are the main sequence progenitors of low-mass AGB stars, show rota- 
tional velocities of a few hundred kms -1 (Royer et al. 2002; Royer, Zorec & Fremat 2003). 
The importance of rotation as a physical process in AGB stars is not restricted to mixing 
and nucleosynthesis. Rotation in stars during their late evolutionary phase possibly drives 
the shaping of bipolar proto-planetary nebulae. According to the interacting-stellar-winds 
model (Kwok 1982; Balick 1987) a fast ionized stellar wind interacts with an equatorially 
dense AGB circumstellar envelope that could be the result of inhomogeneities associated 
with rotation. This model can qualitatively explain the presence of sharp radial structures 
and the wide variety of shapes found in planetary nebulae (see also Icke et al. 1992). Reimers 
et al. (2000) find that elliptical or weakly bipolar planetary nebulae shapes can result from 
dust-driven winds of rotating AGB stars. Soker (2001) considers the possibility that in- 
creased rotation in an AGB envelope may result from swallowing another celestial body, like 
a companion star in a binary system or an orbiting planet. Indeed, this hypothesis gains 
support from the recent discovery of many extra-solar planets and the detection of water 
vapour around evolved AGB stars, possibly due to the presence of comets (Melnick et al. 
2001). Garcia-Segura et al. (1999) present hydrodynamical models in which an equatorial 
density enhancement originates around an intermediate mass single AGB star with a fast 
rotating core that can spin up the extended envelope during mixing events associated with 
the He-shell flashes. 

Langer et al. (1999) evolved a 3M stellar model from the main sequence to the AGB 
phase including the effects of rotation on the stellar structure and mixing. They found 
that rotationally induced mixing at the core-envelope interface after a thermal pulse could 
be responsible for the formation of the partial mixing zone that hosts a 13 C pocket and 
subsequent s-process nucleosynthesis. They also found that mixing in the partial mixing 
layer continues throughout the entire interpulse phase. In this paper we investigate mixing 
and the s-process at the core-envelope interface of this stellar model with rotation. We will 
also present a comparison with models including mixing due to hydrodynamic overshooting 
as well as parametric models that further illustrate the effect of slow mixing of the s-process 
layer during the interpulse period. 

In the current model the s-process occurs in every interpulse-pulse cycle from when 
third dredge-up starts until the end of the AGB evolution. Several stellar parameters that 
are important for the computation of the s-process, e.g. the mass of the intershell, the 
temperature at the base of the convective shell, the overlapping factor between subsequent 
pulses and the third dredge-up, are different at each interpulse-pulse cycle. Detailed s-process 
calculations, such as those of Gallino et al. (1998), take into account these effects. However, 
the features of the 13 C pocket are kept the same in each interpulse, and the changes in the 
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temperature are not large enough to affect the modality of the burning of 13 C in different 
interpulses. The only effect is that in detailed calculations the neutron exposure in the pocket 
slightly decreases with the interpulse number because the amount of s-process material 
increases in the intershell. We perform simulations of the s-process over only one interpulse 
period, which in first approximation represents all the interpulse periods. In § 2 we derive 
the basic properties of the partial mixing zone and interpulse s-process layer from average 
observational features and simplified AGB evolution properties. We describe the stellar 
models and the nucleosynthesis code in § 3. The following section is devoted to our scheme 
for the heavy s-process neutron sink (§4). The properties and the effects of mixing induced 
by hydrodynamic overshoot and by rotation are presented in § 5 and § 6 respectively. Mixing 
for the s-process is further explored with synthetic models in § 7. The particular problem of 
the production of 19 F is addressed in § 8 and we present a final discussion in § 9. 



2. Constraints on the Properties of the Partial Mixing Zone 

The properties of a partial mixing zone that reproduces the s-process features observed 
in AGB stars can be studied in detail with models including the effect of many consecu- 
tive thermal pulses and neutron exposure events (Gallino et al. 1998; Goriely & Mowlavi 
2000). Here we derive some basic constraints on the properties of the s-process zone by 
following a much simpler approach. We consider two s-process indicators: the index s/s Q is 
the overproduction factor of s-process elements with respect to the initial solar value. We 
have used for this index the average of the production factors of Y and Nd. The index 
[hs/ls] = [hs/Fe] - [ls/Fe] monitors the distribution of the s-process elements. We have used 
[ls/Fe]=|([Y/Fe]+[Zr/Fe]) and [hs/Fe] = |([Ba/Fe] + [La/Fe] + [Ce/Fe] + [Nd/Fe] + [Sm/Fe]), where 
square brackets indicate the logarithmic ratio with respect to the solar ratio. Observation- 
ally, the spectroscopic studies of the s-process abundances in evolved low-mass stars of solar 
metallicity can be summarized as < log(s obs /s ) < I and —0.5 < [hs/ls] < (Busso et al. 
1995). 

The observed overproduction factors in the envelope are related to the overproduction 
factors in the s-process zone by two dilution factors that result from two subsequent mixing 
events: the He-flash convective mixing and the third dredge-up. Assuming that no significant 
amount of s-process material is available in the envelope initially and considering only the 
contribution of the s-process in the interpulse, the abundance of any species in the envelope 
after third dredge-up events in m identical TP cycles is related to the abundance in the 
s-process zone (partial mixing zone, PM) by 



Fenv = q m Yi 



PM 



M IS M C 
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where M DUP , Mis and M env are the masses of the dredged-up layer, the intershell zone covered 
by the He-flash convection and the envelope respectively. In low mass AGB stars with core 
masses of ~ 0.6 M these quantities are of the order M DUP ~ 3 • 1O~ 3 M , M IS ~ 1O _2 M 
and M cnv ~ 0.5 M Q . These masses vary from pulse to pulse and are dependent on the core 
mass and on the treatment of mixing. For example in the model sequence with rotation 
Mis = 1 • 10~ 2 M Q at a core mass of M c = 0.746 M Q (Langer et al. 1999) while in the 
sequence with overshooting at all convective boundaries Mi S = 2.4 • 10" 2 M at a core mass 
of M c = 0.628 M (Herwig 2000). The mass M PM refers to the layer of the partial mixing 
zone that at the end of the interpulse phase contains s-process material. In models without 
mixing during the interpulse phase this corresponds to the region of the partial mixing zone 
were initially —3 < logX(H) < —2 (Goriely & Mowlavi 2000). In models with mixing during 
the interpulse like in the case with rotation the extent of the partial mixing zone can only 
be determined at the end of the interpulse when the calculation has yielded the s-process 
nucleosynthesis and mixing result. The factor q describes the effect of overlapping He-flash 
convection zones in subsequent TP (see below). 

Without resorting to the detailed results of full stellar evolution calculations we estimate 
the number m of thermal pulses with dredge-up events that enrich the envelope in a semi- 
empirical way. As we have discussed in the introduction carbon stars are believed to be 
the result of recurrent third dredge-up events in low mass stars with initial ZAMS masses 
predominantly in the range 1.5 < M/ M < 3 (see also Groenewegen et al. 1997). These 
stars end their lives as white dwarfs with masses in the range 0.57 < Mj M < 0.68, 
according to the revised stellar initial-final mass relation of Weidemann (2000). In fact, 
the mass distribution of white dwarfs peaks at or just below 0.60 M (Koester et al. 1979; 
Weidemann & Koester 1984; Bergeron et al. 1992; Napiwotzki et al. 1999). This mass 
distribution is very similar to that of central stars of planetary nebulae (Stasinska et al. 
1997), which are in the evolution phase between AGB stars and white dwarfs. This means 
that the majority of carbon stars must have achieved the necessary abundance enrichment 
through the third dredge-up before or when reaching a core mass of ~ 0.6 M . According 
to the synthetic AGB models of Marigo et al. (1996) significant dredge-up must commence 
at core masses of 0.58 M or even below in some cases, in order to reproduce the carbon 
star luminosity function (Marigo et al. 1999). Therefore the relevant chemical enrichment of 
AGB stars typically occurs within an effective core mass growth of about AM cg = 0.02 M , 
and maybe up to 0.06 M in some cases. For low mass TP-AGB stars the core mass growth 
per TP is about AMh = 6 • 10~ 3 M and with a dredge-up parameter of A = 0.5 about m=7 
(AM cg = 0.02) and possibly up to 20 thermal pulses with dredge-up mixing (if dredge-up 
starts at M c = 0.54 M , AM cg = 0.06) can be considered to be responsible for the abundance 
enrichment. 
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Expression Eq. (1) does not take into account that processed heavy elements accumulate 
in the intershell from one TP to the next because the He-flash convection zone is partly 
overlapping (overlap factor r) with layers that have been swept by the previous He-flash 
convection zone. In the case of nucleosynthesis of a species (s) during the interpulse phase in 
a partial mixing zone the production stays roughly constant from TP to TP. The abundance 
due to such a production is X s = XpmMpm/Mis- The total intershell abundance at the n th 
TP with third dredge-up is then given by X n = X s + rX n _ x . The number I of TPs needed 
to approach a 90% level of some asymptotic value for the intershell abundance of species s 
is given by I = — 1/ logr. For n > I it follows that 

i-i 

X n ^X s -q = X s -Y j r t (2) 

i=0 

The overlap factor decreases as a function of TP number and also depends on the details 
of the third dredge-up efficiency Typically overlap factors in stellar models decrease from 
about 0.8 at the earliest TPs to an asymptotic value larger than 0.4. The condition n > I 
is approximately satisfied for r ;$ 0.6, for which Eq. (2) returns q = 2.3. Other triplets of 
(r,l,q) are (0.4,3,1.6) and (0.8,11,4.6). By evaluating Eq. (1) for the numbers specified 
above, and with Y cnv in Eq. (1) given by using the maximum log(s bs/s©) = 1, we derive a 
logarithmic expression that relates the s-process overabundance in the PM zone with the 
mass of that zone: 

log M PM = -log(s PM / s Q ) + c (3) 
where c = -0.14 (0, -0.44) for m = 10 (7, 20). 

In Figure 2 we show the variation of log(spM/s©) and [1is/1s]pm with the neutron expo- 
sure. These trends have been obtained by fully implicit network calculations containing the 
s-process nucleosynthesis with neutron captures on all isotopes from He to Pb (as in Lugaro 
et al. 2003b), and all relevant charged particle reactions (see § 3 for reaction rate references). 
As more neutrons are released the average s-process overabundance increases. In the current 
s-process model the partial mixing zone by definition does not extend into the He-burning 
shell, and therefore the partial mixing zone can not exceed the mass of the intershell layer. 
Thus, the mass available for the pocket is Mpm < 10~ 2 M©. This, together with Eq. (3) 
and m — 10 requires that log(spM/s©) > 1-86 and translates into a minimum for r in the 
partial mixing zone. In Figure 2 the left shaded part of the diagram with r < 0.2mbarn _1 
is thereby excluded as the predominant region of r in the partial mixing zone. Values of 
log(spM/s ) < 1.86 in that region would require a prohibitively large partial mixing zone 
(Mp M > 10~ 2 M ) in order to reproduce the most s-process enriched stars of solar metallicity 
with log(s bs/s©) = 1- From [hs/ls] bs < the r- values corresponding to the right shaded 
area in Figure 2 can be excluded as typical for the s-process. In fact, if the neutron exposure 
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values predominant in the pocket exceed ~ 0.5 mbarn^ 1 the envelope s-process abundance 
will show [hs/ls] > 0. This case is discussed in §5.3 with an example in Table 1 (see also 
Lugaro et al. 2003b). To summarize, we estimate the following properties for the partial 
mixing zone in stars of solar metallicity: 

• the 13 C pocket should generate a predominant neutron exposure in the range ~ 0.2 - 
0.5 mbarn -1 ; 

• using Eq. (3) and log(spM/s©) = 4, corresponding to the upper limit of r, the mass 
of the partial mixing zone should obey Mpm > 7 • 10~ 5 M (assuming m = 10), in 
agreement with the detailed studies. 

3. Stellar Evolution and Nucleosynthesis Models 

Stellar evolution models with rotation are computed in the same way as in Langer 
et al. (1999). The lD-hydrodynamical stellar evolution code (Langer 1998; Heger et al. 
2000) considers angular momentum, the effect of centrifugal force on the stellar structure 
and the following rotationally induced transport mechanisms for angular momentum and 
chemical species: Eddington-Sweet circulation, Solberg-H0iland and Goldreich-Schubert- 
Fricke instability, and dynamical and secular shear instability. The /i-gradient, which acts as 
a barrier for rotationally induced mixing, and the Ledoux criterion for convection and semi- 
convection are considered. The nuclear energy generation is computed in the operator split 
approximation with a network including pp-chain, H-burning cycles and He-burning. For 
comparison we analyze stellar models with hydrodynamic overshoot of Herwig (2000). These 
models feature a partial mixing zone of protons and 12 C very similar to that assumed for 
s-process calculations with a parameterized partial mixing zone, like in Goriely & Mowlavi 
(2000). They assume however that overshoot is to some extent present at all convective 
boundaries. The two processes of rotation and overshoot are considered independently from 
each other as the stellar evolution models with rotation do not include overshoot and the 
overshoot models do not include rotation. 

The nuclear reactions under consideration for the s-process do not contribute in a signifi- 
cant way to the energy generation in the star. Hence post-processing is a valid approximation 
and a faster and easier-to-handle alternative to re-computing the whole stellar evolution in- 
cluding all the needed nuclear species. An important feature of our rotating AGB stellar 
models is a weak and persistent mixing of the partial mixing layer where the 13 C and 14 N 
pocket are located next to each other. The s-process can not be computed under the as- 
sumption of stratification anymore, as done for example by Lugaro et al. (2003b). Therefore 
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we have developed a nucleosynthesis and mixing post-processing code (SBM6) which solves 
simultaneously for mixing processes, charged particle nuclear reactions as well as neutron 
production and destruction reactions and /3-decays. The code uses thermodynamic and mix- 
ing properties of the stellar evolution models as input and follows the abundance evolution 
of the nuclear network described below. In the stellar evolution code with rotation the time 
and spatial resolution is determined by the H-burning shell. Grid rezoning after the ther- 
mal pulse at the core-envelope interface may introduce some numerical diffusion which can 
be avoided if the s-process simulation is carried out on a fixed grid beginning right after 
the formation of the partial mixing zone. In the post-processing we use a high-resolution, 
equidistant Lagrangian grid with 400 points to cover our partial mixing layer which has a 
mass of 1O~ 4 M . The solution scheme of coupled burning and mixing is fully implicit (Her- 
wig 2001) with adaptive time steps. A Newton-Raphson iteration is accepted as solution if 
the corresponding greatest relative correction is less than 10~ 3 for any species at any grid 
point with a mass fraction X > 10~ 21 . This guarantees a numerical precision of 0.1% for 
any species at each time step and the correct integration of the neutron abundance, which 
is treated like that of any other species. The neutron density typically encountered during 
the interpulse (N n < 10 7 n/cm~ 3 ) corresponds to X(n) < 2 • 10~ 21 . In practise the precision 
for neutrons is much better than 0.1% even for X(n) < 10~ 21 because the convergence of 
the set of equation is determined by other species. Due to the fully implicit solution scheme 
only a fraction (about one hundredth) of the time steps of the evolutionary code need to 
be computed with the SBM6 code because in the partial mixing zone the thermodynamic 
conditions change much slower than at the location of the H-burning shell. 

For the post-processing simulations we have considered the following reactions in the 
SBM6 code: 

• (p,7) from NACRE compilation (Angulo et al. 1999) on 12 C, 13 C, 13 N, 14 N, 15 N, 16 0, 
17 0, 18 0, 19 F, 20 Ne, 21 Ne, 22 Ne, 22 Na, 23 Na, 24 Mg, 25 Mg(both to 26 Al g and 26 Al m ), 
26 Mg, 27 Al, 28 Si, 29 Si and 30 Si. 31 P and 26 Al g (p, 7 ) from Iliadis et al. (2001). 14 C(p, 7 ) 
from Wiescher et al. (1990). 

• (p, a) from NACRE on 15 N, 17 0, 18 and 19 F. 

• /3-decays have been assumed to follow instantaneously were applicable, except for the 
/5-decays of 22 Na, 26 Al g , 13 N, 14 C and 59 Ni from the Karlsruher Chart of Nuclides. 

• (a, 7) from NACRE on 12 C, 15 N and 18 0, and on 14 C from Gai et al. (1987); Hashimoto 
et al. (1986); Funck & Langanke (1989) 1 . 



We have used the electronic NETGEN database to retrieve these and other rates (Jorissen & Goriely 
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• (a, n) from NACRE on 13 C and 22 Ne. 

• (n, 7) from the compilation by Lugaro (2001, Appendix B) which is largely based on 
Bao et al. (2000) on 12 C, 14 N, 16 0, 21 Ne, 20 Ne, 22 Ne, 23 Na, 24 Mg, 25 Mg, 26 Mg, 27 Al, 
28 Si, 29 Si, 30 Si, 31 P, 56 Fe, 57 Fe, 58 Fe, 59 Co, 58 Ni, 59 Ni, 60 Ni, 61 Ni and 62 Ni. 

• (n, p) neutron recycling reactions on 14 N (Gledenov et al. 1995) and 59 Ni from the 
Reaclib Data Tables of nuclear reaction rates (1991), updated version of the compilation 
by Thielemann et al. (1987). 

• A light and a heavy neutron sink accounts for neutron absorbing species between 32 S 
and 56 Fe and above 62 Ni respectively (see §4 for details). 

We consider the light n-capture reactions that are important to determine the n-density. 
The efficiency of isotopes to absorb neutrons can be measured by the product of the n- 
capture Maxwellian-averaged cross section and the number density of the isotope. Assuming 
a solar abundance distribution the 12 most important light n-absorbing isotopes are (in 
the order of decreasing absorbing efficiency, (n, 7) reaction unless noted): 14 N(n,p), 30 Si, 
27 Al, 28 Si, 29 Si, 25 Mg, 16 0, 14 N, 24 Mg, 23 Na, 20 Ne and 12 C. The order of this list of n- 
capture isotopes changes at locations in stars where a non-solar abundance distribution has 
been established. E.g. in the 14 N pocket (see § 5) the 14 N abundance is so high that the 
n-captures by 14 N(n, 7) reactions outnumber all other (n, 7) reactions on light elements and 
is itself dwarfed by 14 N(n, p) 14 C. The long-lived ground state of the unstable nucleus 26 Al is 
also a major neutron poison due to its very large cross sections for (n, p) and (n, a) reactions 
(Skelton et al. 1987; Koehler et al. 1997). This nucleus is produced by proton captures on 
25 Mg during H-burning, however it is efficiently destroyed by neutron captures in the thermal 
pulse convective regions, hence it is not present in the 13 C pocket and we have omitted it in 
the computation. 

The neutron recycling induced by the 14 N(n,p) 14 C reaction is of particular importance 
for the s-process in rotating AGB stars. This reaction deprives the Fe-seed and trans-iron 
elements of available neutrons. The major effect is that the neutron exposure is smaller if 14 N 
is present (see Appendix A). In addition the 14 N(n,p) 14 C reaction opens a channel to make the 
isotope 15 N, hence potentially 19 F. 14 C has a half life of 5730 yr. Under the typical conditions 
of n-production in the interpulse phase (T = 9 • 10 7 K,p = 3700gcm~ 3 ) 14 C can capture an 
a particle or a proton. The 14 C(o;,7) 18 reaction is 2.5 times slower than 13 C(o;, n) 16 0, 
but still affected by large uncertainties. The p-capture reaction of 14 C ( 14 C(p, 7) 15 N) is 3.7 
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times faster than 12 C(p,7) 13 N and about as fast as 13 C(p,7) 14 N. No other p- or a-capture 
reactions of 14 C are important and also the 14 C(n, 7) 15 C can be neglected. 

4. Approximating the s-Process with a Neutron Sink 

We introduce two reactions and two artificial particles that together give an estimate on 
the overproduction of s-process elements and the s-process distribution. First, we consider 
the reaction 62 Ni(n, 7) 63 G, where 63 G is an artificial particle with mass number 63. We 
identify the number abundance of 63 G with the combined number abundance of all isotopes 
heavier than 62 Ni: 

209 

F( 63 G) = r ( As )> ( 4 ) 

A=63 

where S stands for the element symbol of the respective species with A > 62. For a solar 
composition the number abundance of 63 G is Y Q ( 63 G) = 5.4- 10~ 8 . Starting on 56 Fe a suitable 
neutron exposure will quickly lead to the formation of heavy particles A S (with A > 62), 
hence the abundance of 63 G increases. In order to count the number of neutron captures 
that occur on species with A > 62 we introduce a second reaction 63 G(n, 1 L) 63 G, which plays 
the role of heavy neutron sink. The Maxwellian-averaged cross section of 63 G is computed 
in the usual way (Jorissen & Arnould 1989): 

209 

o-( 63 G) = Y^G)- 1 *aY( a S). (5) 

A=63 

Neutron captures will occur repeatedly on individual 63 G particles, thereby simulating the 
production of increasingly heavy s-process isotopes. These neutron captures are responsible 
for the final s-process element distribution of the species represented by 63 G. They do, 
however, not change the number abundance of 63 G. The ratio (L/G) = y( 1 L)/y( 63 G) is 
similar to the customary quantity n cap (Clayton 1968) and is a measure of the s-process 
distribution (Figure 3). In fact, if one defines the artificial sink particle G to be the product 
of a n-capture on 56 Fe, then (L/G) = n cap . With our choice of G being the product of a n- 
capture on 62 Ni, n cap is slightly larger than (L/G) for a given neutron exposure because n cap 
takes into account the n-captures on isotopes from 56 Fe to 62 Ni. These n-captures mainly 
take place during the initial phase of the neutron exposure phase for r < O.lmbarn -1 . 

The neutron cross section of the sink particle 63 G depends on the abundance distribution 
of the particles it represents. In Figure 4 the Maxwellian-averaged sink cross section at 8keV 
according to Eq. (5) is shown during calculations of the s-process nucleosynthesis with the 
network including neutron captures on all isotopes up to Pb, starting with a solar abundance 
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distribution of trans- iron elements. The variation of <r( 63 G) reflects the changing abundance 
distribution of heavy elements. Initially, species accumulate at the stable neutron magic 
nuclei on the s-process path ( 88 Sr, 138 Ba, 208 Pb) and the averaged cross section decreases. 
As more neutrons are released the contribution of species in between increases and finally 
dominates the sink cross section. The choice of a constant sink neutron cross section for 
the entire s-process simulation is the largest individual uncertainty when approximating the 
s-process with artificial sink reactions as described above. However, it turns out that the 
error introduced by using a constant sink cross section is sufficiently small for our purpose. 
In Figure 5 we show the neutron exposure of calculations with the sink treatment with three 
choices of the sink cross section as compared with s-process calculations performed with a 
network including neutron captures on all isotopes up to to Pb for the three cases of Figure 4. 
For low neutron exposures the influence of the sink cross section is small, because the neutron 
density is dominated by 56 Fe and the lighter neutron capture elements. For the s-process in 
the partial mixing zone of stars of solar metallicity the most important range for r is between 
~ 0.2 . . . O.Smbarn" 1 . We choose for our simulations er( 63 G) = 120mbarn which reproduces 
the neutron exposure from the s-process calculation up to Pb within 10%. 

We also introduced a light neutron-sink reaction to take into account neutron captures 
on elements from S to Mn, which are missing in our network. For this reaction, which plays 
only a minor role, we have used the value <7ii g ht(8keV) = 7.36mbarn given by Lugaro (2001). 

5. Mixing from hydrodynamic overshoot 

The hydrodynamical properties of convection inevitably result in some turbulent mixing 
into the stable layers adjacent to convectively unstable regions. In fact, any model of con- 
vection in the hydrodynamical framework predicts that the turbulent velocity field decays 
roughly exponentially inside the stable layers (e.g. Xiong 1985; Freytag et al. 1996; Asida 
& Arnett 2000). A depth- and time-dependent hydrodynamic overshoot approximation has 
been used in a number of recent studies (Herwig et al. 1997; Mowlavi 1999; Mazzitelli et al. 
1999; Salasnich et al. 1999; Herwig 2000; Cristallo et al. 2001), with the aim to capture the 
main consequences of hydrodynamic mixing into the stable layers induced by convection. 

5.1. Third dredge-up and H/ 12 C-partial mixing zone 

It has been shown by Herwig (2000) that overshoot at all convective boundaries, in- 
cluding the base of the convective envelope and the bottom of the He-flash convection zone, 
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strongly increases the efficiency of the third dredge-up at low core masses. This is required 
observationally in order to reproduce the observed C-star luminosity function in the Magel- 
lanic Clouds (Richer 1981; Frogel et al. 1990). Synthetic models of the AGB phase in which 
the third dredge-up parameter is derived observationally have demonstrated that typically 
efficient third dredge-up must take place at core masses as low as 0.58 M (Marigo et al. 
1996) or at even lower core masses around 0.54 M Q according to a more recent analyzes by 
Marigo et al. (1999). This condition is not met by most stellar models, including those of 
Mowlavi (1999) who considers hydrodynamic overshoot only at the bottom of the convective 
envelope (see his Fig. 10b for a comparison of third dredge-up efficiencies found by different 
authors). The model grid of Herwig et al. (2000), which includes overshoot at all convective 
boundaries, includes cases (e.g. the 2M , Z=0.01 case) which cover the low C-star lumi- 
nosity tail as required by observations. The third dredge-up properties of AGB models are 
important for the s-process because third dredge-up is needed to bring the processed material 
to the surface. Hydrodynamic overshoot creates a partial mixing zone at the core-envelope 
interface with a continuous decrease of the H/ 12 C-ratio from the envelope into the intershell 
layers (Herwig et al. 1997; Herwig 2000). 

The global properties of the s-process can be reproduced with a partial mixing zone 
resulting from hydrodynamic overshoot (Goriely & Mowlavi 2000; Lugaro & Herwig 2001). 
The main features of the s-process overabundance distribution are mainly determined by the 
regions of the pocket which have the largest neutron exposures, not so much by the detailed 
shape of the H/ 12 C-profile within the partial mixing zone. Even if the treatment of hydro- 
dynamic overshoot according to an exponentially decaying velocity field is not correct and 
the actual functional dependence of overshoot efficiency with depth is somewhat different, 
then the s-process will most likely be affected only slightly as long as the H/ 12 C-profile is 
somehow continuous. For example, Denissenkov & Tout (2002) investigate gravity waves be- 
low the convective envelope as a cause for extra mixing to produce partial mixing of protons 
and 12 C. This mixing process is another way of looking at the mixing resulting from the 
perturbation of the convective boundary due to turbulence and leads to neutron exposures 
in the region close to that of previous models featuring a continuous decrease of the proton 
abundances into the 12 C-rich core. 

An additional effect is introduced in models that consider overshoot at all convective 
boundaries (Herwig 2000). In these models the 12 C abundance in the intershell is about 
twice as large than in models without overshoot at the base of the He-shell flash convection 
zone. Lugaro et al. (2003b) have shown that the neutron exposure in the s-process layer is 
proportional to the 12 C abundance in the intershell. Hence, in models that consider overshoot 
at the base of the He-shell flash convection zone the neutron exposure in the s-process layer 
is higher than in models that do not include this overshoot. 
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5.2. How much overshoot? 

The initial computations of AGB stars with hydrodynamic overshoot were carried out 
with an exponential overshoot parameter of / = 0.016, which was motivated by the efficiency 
derived from convective core overshoot of main sequence stars. The effective mass of the 
partial mixing zone where the neutrons are efficiently released is confined within the region 
where the proton abundance follows —2 < logX P < —3 (for an intershell 12 C mass fraction 
of ~ 20%, Goriely & Mowlavi 2000). According to this criterion the mass of the s-process 
layer computed with / = 0.016 is only ~ 1O _6 M , which is much smaller than required (see 
§2). 

However, one overshoot efficiency parameter will not apply to all convective bound- 
aries during all evolutionary phases. After Shaviv & Salpeter (1973) first considered the 
possibility of convective overshoot several studies have used a very simple prescription in 
which convective mixing was treated instantaneously and overshoot was simply a matter of 
extending the instantaneously mixed region by some fraction of the pressure scale height. 
In this approximation main sequence core overshoot should extend by about 0.2H p (e.g. 
Schaller et al. 1992, and references there). Alongi et al. (1991) argued that overshoot of 
0.7H P below the envelope of red giant stars could align the location of luminosity bump with 
observations. The 2D radiation-hydrodynamic simulations by Freytag et al. (1996) have 
shown that the shallow surface convection zone of white dwarfs have exponential overshoot 
mixing according to an overshoot parameter of / = 1.0 while the convection zone simulation 
of A-stars show / = 0.25. For the oxygen burning layer in pre-supernova models Asida & 
Arnett (2000) found perturbations of the stable layers reaching one pressure scale height 
beyond the formal convective boundary. Thus, there is ample indication that the overshoot 
efficiency is not the same at different convective boundaries. However, convective overshoot 
is not a stochastic process as long as the convective turn-over time scale is shorter than the 
thermal time scale of the region which hosts the convective boundary. For similar convective 
boundaries one should expect a similar overshoot efficiency. This expectation is supported 
by 2D hydrodynamical computations by Deupree (2000) who showed that the core overshoot 
distance of ZAMS stars varies only mildly with stellar mass. 

Here we choose an exponential overshoot parameter for the hydrodynamic overshoot at 
the bottom of convective envelope of / = 0.16. This larger overshoot is only applied during 
the third dredge-up phase. This has no major effect on the properties of the models, other 
than stretching the partial mixing zone and consequently the 13 C and 14 N profiles in that 
layer over a larger mass range. The peak neutron exposure and the s-process abundance 
distribution in the partial mixing layer are not much changed. As a side effect the third 
dredge-up efficiency is slightly increased, by ~ 20%. 
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As mentioned in the introduction observations as well as the analysis of pre-solar mete- 
oritic SiC grains suggest that stars with otherwise identical initial conditions have a range of 
s-process efficiencies. Such a range can not be expected to result from overshoot since such 
mechanism is not expected to be a stochastic process. 

5.3. Neutron production for the s-process in the overshoot model 

We model the abundance evolution in the partial mixing zone during the 7 th interpulse 
phase of the 3 M , Z = 0.02 sequence of Herwig et al. (2000), with an overshoot efficiency / = 
0.16 during the third dredge-up phase. We use the post-processing code (SBM6) described 
in §3 because the computation of the stellar evolution does not include all the species and 
reactions needed to study the s-process. As initial conditions we use the thermodynamic 
and abundance profiles from the stellar evolution model at the end of the third dredge-up 
phase after the thermal pulse. These profiles are mapped to the equidistant, Lagrangian 
post-processing grid and then evolved according to the stellar structure models at a series 
of times throughout the interpulse phase. 

We start the simulation with the partially mixed H/ 12 C zone of ~ 10~ 4 M that has 
formed at the end of the third dredge-up phase as a result of time- and depth-dependent 
hydrodynamic overshoot (top panel, Figure 6). In this model no mixing takes place during 
the interpulse phase. In the middle panel of Figure 6 the 13 C neutron source has started 
releasing neutrons and up to 10% of the 13 C abundance has been consumed. In the upper 
part of the partial mixing zone where 14 N dominates the majority of neutrons is absorbed 
by the 14 N(n,p) 14 C reaction. This can be seen from the profile of 14 C. The maximum 
neutron density is 5.7 • 10 6 cm -3 at that time and it subsequently reaches a peak value of 
1.1 • 10 7 cm~ 3 . The bottom panel of Figure 6 shows the total neutron exposure in the s-process 
layer before it is engulfed by the subsequent He-flash convection zone. A maximum neutron 
exposure r = 0.72mbarn~ 1 is found. This is about a factor two larger than in previous 
s-process simulations of the partial mixing scenario (Goriely & Mowlavi 2000) because in 
our overshoot models the 12 C intershell abundance is larger than in models that do not use 
overshoot at the bottom of the He-flash convection zone. 

In comparison to the observational constraints for the partial mixing and s-process 
layer derived in § 2 we find that the neutron exposure in this model with hydrodynamic 
overshoot at all convective boundaries is too large. Lugaro et al. (2003b) concluded that 
to match observed stellar s-process abundance patterns the overshoot at the bottom of the 
He-flash convection zone should be smaller than that used by Herwig (2000, f=0.016). For 
example, the higher temperatures in the He-flash convection zone due to overshoot result 
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in very large 96 Zr/ 94 Zr ratios, in contradiction with measurements in mainstream pre-solar 
meteoritic SiC grains. The temperature in the He-flash convection zone decreases with the 
stellar mass. However, some SiC show almost no presence of 96 Zr and they are difficult to 
explain by current models even by considering stars with mass as low as 1.5 M (Lugaro 
et al. 2003a) and no overshoot included. The inclusion of overshoot at the base of the 
He-flash convection zone increases the temperature for a given stellar mass, making it even 
more difficult to explain these data. On the other hand, current models of hydrogen-deficient 
post-AGB stars (Herwig et al. 1999) can only reproduce the abundance patterns observed in 
PG1159-type objects and central stars of planetary nebulae of spectral type [WC] (Koesterke 
& Hamann 1997; Dreizler et al. 1996) with AGB progenitor models including intershell 
overshoot. Further investigations of extra mixing have to address whether the constraints 
from s-process branchings and from subsequent evolutionary stages can be resolved in well 
adjusted models. The mass of the layer in which significant overproduction factors of 
s-process elements can be expected is in our calculation of the order ~ 3 • 10 -5 M Q , about a 
factor of two less than the minimum pocket mass required according to the estimates in §2. 
The choice of the overshoot parameter at the base of the convective envelope only affects 
the extent in mass of the partial mixing zone, and very little the s-process distribution. The 
latter is determined by the neutron exposure in the pocket. For that reason one can not 
remove the problem of the large neutron exposure in our overshoot model by a fine tuning 
of the overshoot parameter at the base of the convective envelope. 

This can easily be shown with a numerical test carried out with a parametric s-process 
code, as that used by Lugaro et al. (2003b, and ref. there). We made four runs of the 
s-process model with radiative 13 C burning through 13 thermal pulse events and mixing 
episodes according to the stellar evolution model including overshoot, as described in Lugaro 
et al. (2003b). To select the effect of varying of the 13 C-pocket mass only, we neglect the 
contribution of the 22 Ne neutron source during the He-shell flash. The benchmark (BM) 
case corresponds to the 13 C pocket obtained with an exponential overshoot parameter of 
/ = 0.128 at the bottom of the convective envelope. The pocket is computed starting from 
an intershell 12 C mass fraction of 0.43 and is kept identical in all the interpulse phases. We 
multiply/divide the mass of the pocket as indicated in Table 1, which approximates the effect 
of larger /smaller / values. Note that this is different than changing the 13 C abundance in 
the pocket as done in order to obtain a spread in neutron exposures by Busso et al. (2001). 
Table 1 gives s-process indices at the stellar surface at the end of an AGB model sequence. 
In all cases [hs/ls]> 0, clearly in contrast to the observed negative [hs/ls]. The smallest 
[hs/ls] values which are closest to the observed values correspond to almost no s-process 
overproduction ([ls/Fe] and [hs/Fe] ~ 0) and must be excluded. 
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6. Mixing in rotating AGB stars 

For the formation of the 13 C pocket and the s-process in TP- AGB stars the most relevant 
rotationally induced mixing instability is caused by shear at locations of large angular velocity 
gradients. In order to anticipate rotationally induced mixing events it is therefore important 
to consider the angular velocity evolution throughout the thermal pulse cycle. Throughout 
this and the following section we make use of the interpulse phase = t -j^ L , where t is 
the time of the thermal pulse, i.e. the peak He-burning luminosity and At is the interpulse 
period of 3 • 10 4 yr. 

The evolution of angular velocity in the intershell region during the 25 th interpulse 
period of our 3 M rotating TP- AGB model star is shown in Figure 7. The main features of 
the angular velocity evolution can be understood by redistribution of angular momentum in 
the convective regions and conservation of angular momentum in contracting or expanding 
stable layers. The first profile at = —0.0536 shows a model of stable H-burning during 
the interpulse phase shortly before the He-shell flash. The H-shell is located just below 
the mass coordinate 0.746 M Q . The steep jump at ~ 0.741 M marks the largest extent of 
the He-flash convection zone during the previous thermal pulse. The He-flash convection 
zone redistributes angular momentum from the faster rotating C/O-core throughout the 
intershell and the layers just below the H-rich envelope. The profile shows a plateau between 
0.741 and 0.746 M Q . This region contains the H-shell ashes which have been deposited there 
by the outward burning H-shell. The angular velocity is larger in this area than in the 
convective envelope above the H-burning shell. Angular momentum is efficiently distributed 
by convection throughout the convective envelope, which rotates almost rigidly. However, 
the H-burning ashes contract onto the core and therefore the angular velocity increases, 
hence the plateau between 0.741 and 0.746 M is formed. 

At = 0.0034 the He-flash convection zone has formed and the mass layers covered 
by convection are spun up. At this time the He-flash convection zone reaches up to mass 
coordinate 0.745 M . Between this and the following profile at = 0.0037 the He-flash 
convection zone reaches its fullest extent up to 0.746 M (just below the location of the 
H-burning shell before the thermal pulse). In addition, the intershell region is expanding 
due to the energy provided by the He-flash. This leads to a reduction of the angular velocity 
in the intershell. This effect is still present at = 0.0087, until contraction resumes and the 
angular velocity in the intershell region increases accordingly (profiles = 0.0122 to 0.1706). 
The angular velocity profiles at later times show again the formation of a plateau between 
0.746 and 0.749 M , where the H-shell is now located, due to deposition of H-shell ashes on 
top of the core. 

From the evolution of the angular velocity in rotating TP-AGB stars we can anticipate 
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the following mixing properties. Efficient shear mixing will take place at ~ 0.746 M after the 
formation of the large angular velocity gradient. This gradient forms because two convective 
regions extend to this mass coordinates from below and above in short succession. The 
He-flash convection first taps the reservoir of high angular momentum in the core. Then, 
the convective envelope establishes contact between the fast-rotating intershell and the low- 
rotating envelope. At this interface rotationally induced shear mixing will be most efficient in 
acting on abundance gradients. The timescale to establish a partial mixing zone of hydrogen 
and 12 C as needed for the formation of the 13 C pocket is limited to the time interval when 
the envelope and the intershell are in contact, i.e. between the end of the third dredge-up 
and the reignition of the H-shell. This period lasts 2000 to 3000 yr (Figure 1). However, as 
can be seen in Figure 7 the steep angular velocity gradient remains at this mass coordinate 
after the formation of the partial mixing zone and shear mixing at this location will persist 
throughout the intershell period. 

Correspondingly, two mixing periods during the interpulse phase resulting from rotation 
can be distinguished as shown by the temporal evolution of the mixing coefficient presented 
in Figure 8. During the initial envelope-core contact period at the end of the third dredge-up 
a partial mixing zone of protons and 12 C forms (top panel, Figure 9). The mixing efficiency 
in Eulerian coordinates at this time is log A- ~ 5 (in cgs units), but decreasing rapidly at a 
time scale of ~ 1000 yr. However, after the thermal pulse the mass gradient in the partial 
mixing zone increases steadily as the intershell layers contract and evolve towards a pre- 
flash structure. For assessing the effect on the abundance evolution the Lagrangian mixing 
coefficient should be considered: 



where all symbols have their usual meaning. 

In Table 2 we give the mixing properties in the partial mixing zone at three different 
times: an early phase soon after the pulse, a phase just before the release of neutrons 
starts, and a phase after the neutron source 13 C is exhausted. Despite the comparatively 
large mixing coefficients rotationally induced mixing is not important immediately after the 
initial formation of the partial mixing zone. By the time the neutron source is about to be 
activated the mass gradient has increased by more than two orders of magnitude, which more 
than offsets the decrease of the geometrical mixing coefficient. The Lagrangian mixing range 
then exceeds the mass of the initial partial mixing zone. This trend continues throughout 
the time of occurrence of the s-process during the interpulse phase, as the intershell layer 
and the partial mixing zone gradually contract. 

Finally, we note that the present model with rotation does not show sufficient third 
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dredge-up. This has already been reported by Langer et al. (1999) and is evident from their 
Fig. 3 and 4. We found that the third dredge-up in models of rotating AGB stars is related to 
the efficiency of /^-gradients to inhibit mixing. This is a somewhat uncertain parameterized 
property of rotating stellar models and a different approach in describing this effect could 
change the resulting third dredge-up. 

6.1. Neutron production for the s-process with rotation 

We model the abundance evolution in the partial mixing zone during the 25 th interpulse 
phase of a stellar model of 3 M , Z = 0.02, and initial equatorial rotation velocity of 250km/s 
(Langer et al. 1999). The core mass at this interpulse phase is M c = 0.746 M and the third 
dredge- up only just dips into the carbon-rich intershell. To compute the s-process we used 
the post-processing code (SBM6) described in §3. As in the computation with overshoot, 
we use as initial conditions the thermodynamic and abundance profiles from the stellar 
evolution model at the end of the third dredge-up. We map these profiles to the equidistant, 
Lagrangian post-processing grid and then evolve them according to the stellar structure 
models at a series of times throughout the interpulse phase. As a result of shear mixing 
due to the large angular velocity gradient immediately after the end of the third dredge-up 
a partial mixing zone forms (top panel, Figure 9) similar to that in models with diffusive 
overshoot (Figure 6). The partial mixing zone defined as having a proton abundance of 
-2 < logX(p) < -3 has a mass of M PM = 6 • 1O~ 6 M . 

As described above, throughout the interpulse period shear mixing at the former core- 
envelope interface continues at a low level. This has two effects on the 13 C initially formed 
in the partial mixing zone. The 13 C is spread out, hence diluted, over the mass range 
subject to shear mixing. Secondly, the 14 N and the 13 C pockets, which are well separated 
in the overshoot model, are effectively mixed. As a result the 14 N abundance exceeds the 
13 C abundance everywhere in the pocket (middle panel, Figure 9). In this situation most 
neutrons released by 13 C(a;,n) are absorbed by the reaction 14 N(n, p) 14 C rather than by heavy 
nuclei. Due to the increased amount of neutron absorbers the neutron density is drastically 
reduced; the maximum neutron density is iV n = 2.6 • 10 6 cm~ 3 . This corresponds to the very 
small neutron exposure obtained in this model (bottom panel, Figure 9). The activation 
of the 14 N(n, p) 14 C reaction can be seen from the 14 C profile in the middle panel. Protons 
released by this reaction are captured by 12 C, building more 13 C, as well as by 13 C itself (and 
by 18 0, see discussion in §8). 

In comparison to the basic requirements for the s-process layer derived in § 2 we conclude 
that this particular model of rotational mixing in the TP-AGB interpulse period can not 
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produce the s-process abundance patterns observed in stars of solar metallicity. The neutron 
exposure reaches only r = 0.04mbarn~ 1 , too small for any significant s-process enhancement. 
As far as the mass of the s-process layer is concerned we note that shear mixing has caused 
a substantial broadening. In the end, a zone of 5 • 10~ 5 M has experienced a low neutron 
irradiation. 

One peculiarity of s-process simulations including mixing throughout the interpulse 
period concerns the sensitivity to uncertainties of the nuclear reaction rates. In the param- 
eterized s-process without mixing during the interpulse each layer is treated individually 
and all 13 C is burnt before the onset of the next He-shell flash. As a result, for example a 
different 13 C(o;,n) 16 rate, as long as it allows all 13 C to be consumed during the interpulse 
period, would lead to a slight time shift of the neutron release, but in the end this does not 
affect the neutron exposure. The case with mixing throughout the interpulse is different due 
to the more important recycling of neutrons by 14 N(n, p) 14 C and 12 C(p, 7) 13 N. For example, 
we find that the neutron exposure increases by about 20% when using the NACRE (Angulo 
et al. 1999) 13 C(a,n) 16 rate instead of the CF88 (Caughlan & Fowler 1988) rate. 

In summary, we have demonstrated the main properties of partial mixing layers for the 
s-process due to rotationally induced mixing and due to hydrodynamic overshoot. Both 
models do not seem to be consistent with the observed properties. The overshoot model 
shows some problems related to the overshoot at the base of the He-flash convection. The 
model including rotational mixing shows that the s-process during the interpulse is easily 
prevented if the delicate process of 13 C formation and neutron release is disturbed by mixing. 
In fact, our preliminary model of an AGB star evolving from a main sequence star rotating 
more slowly than the model presented above indicates that also in this case the s-process 
nucleosynthesis during the interpulse period is very much inhibited by the rotation induced 
mixing. 



7. Synthetic Mixing Models 

In the previous section we have encountered the limitations of both mixing processes, 
hydrodynamic overshoot and rotationally induced mixing. Here we want to explore the 
possibility that some combination of both processes may provide a satisfactory mixing law 
for the s-process. To that end we construct some synthetic post-processing models that 
include overshoot and rotation in a simple parametric scheme. 

We approximate the details of rotationally induced mixing by assuming a constant 
mixing coefficient Dip throughout the interpulse phase for the post-processing simulation of 
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the s-process layer. To show that this approach leads to the same results as the detailed 
computation presented in §6.1 we compare the evolution of the neutron exposure with time 
from the detailed computation with that obtained in two test cases (Figure 10). Test rotl 
is computed by starting from the partial mixing zone obtained in the rotation sequence (top 
panel, Figure 9), but excluding any further mixing during the interpulse. The maximum 
neutron exposure is r « 0.38mbarn _1 , much higher than the case when mixing continues 
during the interpulse and very similar to that found in other partial mixing models without 
rotation. However, the neutron exposure integrated over the simulated mass range is two 
times smaller in test case rotl compared to the case when the further mixing during the 
interpulse is included. This reflects the higher degree of n-recycling in a region of wider mass 
in the case when rotation mixing during the interpulse is included. Test rot 2 is computed by 
starting again from the partial mixing zone obtained in the rotation sequence and applying a 
constant Eulerian mixing coefficient log Dip = 1 (cgs units) across the simulation range. This 
procedure can approximately reproduce the key quantities of the detailed post-processing 
simulation presented in §6.1: Mmp < 1CT 4 M Q and r max ~ 0.03mbarn _1 . The fact that the 
test case rot2 approximates the results of the detailed simulation confirms that a constant 
mixing coefficient throughout the interpulse period mimics the second phase of rotationally 
induced mixing that is responsible for the admixture of the 14 N neutron poison and the 
dilution of the 13 C pocket. 

With this simple representation we study the effect that mixing of the s-process layers 
during the interpulse might have on the partial mixing zone formed as a result of hydrody- 
namic overshoot. We perform another set of simulations in which we start with the partial 
mixing layer produced from hydrodynamic overshoot (shown in the top panel of Figure 6). 
A constant mixing coefficient is imposed during the interpulse simulation. Cases for four 
different interpulse mixing efficiencies are shown in Figure 11. Even with very weak constant 
mixing of log Dip = —2 the peak neutron exposure is somewhat reduced compared to the 
overshoot case without interpulse mixing shown in Figure 6. For faster interpulse mixing the 
s-process layer becomes broader and the final neutron exposure declines. 

From these tests we tentatively propose a two-step-mixing scheme for the s-process 
in AGB stars. At the end of the third dredge-up a fast mixing process induced by the 
convectively unstable envelope leads to the formation of a partial mixing zone. A mass 
range of « 5 ■ 1CT 5 M should have a proton abundance of —2 < logX(p) < —3. The 
intershell 12 C abundance should be larger than in models that confine mixing in the He- 
shell flash convection zone strictly to within the Schwarzschild boundaries. Due to mixing 
during the interpulse phase the neutron release is modulated (in any case reduced) and 
a set of otherwise identical stars with a variety of mixing efficiencies at the core-envelope 
layers will display a spectrum of neutron exposures. In this scenario involving hydrodynamic 
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overshoot as well as interpulse mixing, the efficiencies needed to cover the spread in s-process 
efficiencies observed in stars (Busso et al. 2001) and pre-solar grains (Lugaro et al. 2003a) 
can be achieved with mixing efficiencies of the order log Dip ~ . . . — 1.3. This range is 
not currently reproduced by stellar evolutionary sequences including rotation, which predict 
much larger mixing efficiencies (see §6.1). Smaller mixing efficiencies would result from a 
smaller angular velocity gradient at the core-envelope interface and/or a smaller rotation 
rate. This could be achieved by a more efficient angular momentum transport during the 
progenitor evolution. Efficient third dredge-up and/or penetration of the convective pulse 
into the core might affect the mixing properties due to rotation as angular momentum could 
be transported out of the core into the envelope and carried away by mass loss. This effect 
is not present in our computation as our models do not show efficient third dredge-up. 
Moreover, in our models we have not considered magnetic fields which could enhance the 
coupling of core and envelope and decelerate the core (Spruit 1998). 



8. Fluorine Production in Models with Rotation 

Jorissen et al. (1992) have observed high enhancements of fluorine in AGB stars. The 
solar abundance is X( 19 F) = 4.1 • 10~ 7 and observationally X( 19 F) = 6.6 ■ 10~ 7 at the start 
of the AGB phase (average of K and M giants). The typical abundance of 19 F in TP-AGB 
stars (S, M and C stars) is X( 19 F) m (4.5 ± 4) • 10~ 6 . In particular the N-type carbon 
stars show a large spread in 19 F abundances within a small range of C/O values. Overall 
thermal pulses appear to cause a tenfold increase in 19 F as a result of thermal pulse and/or 
interpulse nucleosynthesis and third dredge-up. While the observed correlation of fluorine 
enhancement with s-process enhancement is reproduced by current models (see Figure 13 of 
Goriely & Mowlavi 2000), the observed correlation of fluorine with carbon remains unsolved. 

Several nucleosynthesis paths have been considered and the most likely involves a (n,p) 
reaction (Jorissen et al. 1992). Because the 19 F(a, p) 22 Ne rate is about ten times larger 
than the 22 Ne(a, n) 25 Mg rate at T = 3 • 10 s K (Caughlan & Fowler 1988; Angulo et al. 
1999 ) 22 Ne can be excluded as the n-source for the required (n,p)-reactions. Even if fluorine 
is made during the interpulse period, massive AGB stars can be excluded as 19 F producers 
because they have very hot He-flash convective zones, in which 19 F(a, p) 22 Ne is very efficient. 
Moreover, in massive AGB stars, proton captures at the hot base of the convective envelope 
destroy fluorine efficiently. 

The 13 C(o;,n) reaction is activated when the H-burning ashes are engulfed in the convec- 
tion zone. However, the amount of neutrons released by the 13 C from the H-burning ashes 
is not enough to produce the required abundance of fluorine. Another possible site for the 
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production of fluorine is the nucleosynthesis occurring in the partial mixing zone during the 
interpulse periods. However, the partial mixing zone is typically unimportant with regards 
to the production of light elements because its mass is very small. In fact, the inclusion 
of the partial mixing zone in model calculations does not appear to increase the predicted 
surface abundance of 19 F. This conclusion can be drawn when comparing the abundance of 
19 F at a given C/O ratio in Figure 12 of Goriely & Mowlavi (2000), who studied fluorine 
production in AGB stars including a parametric partial mixing zone, and in Figure 14 of 
Mowlavi, Jorissen & Arnould (1996), who did not include a partial mixing zone in their 
study. Here we check whether rotation can improve the match with observation. 

During the interpulse period and in the partial mixing zone about 80% of the 15 N is 
produced via the chain of reactions 

14 N(n,p) 14 C(a, 7 ) 18 0(p,a) 15 N 

and 20% comes from the reaction 14 N(n, 7) 15 N (Figure 12). During the interpulse phase the 
temperature is not high enough for the 15 N(a,7) 19 F reaction and the production of fluorine 
occurs later in the He-flash convection zone (note that the 15 N(a, 7) 19 F rate is about 50 times 
faster than the destruction of 19 F by a-capture in the temperature range of interest). For 
simplicity we assume that all 15 N present in the intershell before the onset of the He-flash 
will be transformed into 19 F and no 19 F will be destroyed so that the intershell abundance 
at the time of the third dredge-up is X IS ( 19 F) oc Mi 57V /M IS , where Mis N = X PM ( 15 N)M PM . 
With Eq. (1), the observed abundances mentioned above, using q=2.3, m=20 and the other 
numbers used in § 2, the total 15 N production in the partial mixing zone during the interpulse 
phase must be of the order Mis N > 1O~ 7 M to cover the observational data. 

In Figure 12 we show some results for the 15 N production in the partial mixing zone. In 
the post-processing model of the sequence including rotation 15 N is produced by the reaction 
channels described above. The absolute amount is much smaller than needed. However, 
compared to a case with the same initial partial mixing zone but no interpulse mixing (test 
case rotl, §7) the production of 15 N is larger by a factor 2... 3. We analyzed the 15 N 
production in the parametric models presented in § 7 that combine an overshooting partial 
mixing zone with interpulse mixing of the s-process layer. We find that the 15 N production 
increases with interpulse mixing reaching 8 • 10~ 8 M Q for Dip = 0. If rotation in AGB stars is 
instrumental in producing 19 F one might expect an anti-correlation of 19 F with the s-process 
index [hs/ls] as well as with the s-process enhancement. Currently available observational 
data do not show this. Renewed observational work on 19 F in AGB stars is needed as well 
as more detailed models. In addition the question of reaction rate uncertainties has to be 
revised, in particular in the case of 14 C(a,7) 18 0. 
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9. Discussion 

We have compared two mixing processes for the s-process in low-mass AGB stars. Mod- 
els that include overshoot at all convective boundaries feature neutron exposures that are 
larger than allowed to reproduce the observed s-process properties. This is due to the fact 
that overshoot at the base of the He-shell flash convection zone leads to intershell dredge-up 
of additional carbon from the core into the intershell layer. Models that do not include 
overshoot at all convective boundaries (and include the partial mixing zone in a parametric 
way) feature neutron exposures that can explain only a small fraction of the observations. 
To explain the whole range of observational data some spread in the efficiency of the neutron 
release is required in the scenario of 13 C(a,n) occurring in the interpulse periods. If rotation 
is included the s-process layer located at the core-envelope interface is continuously mixed 
throughout the interpulse period. The s-process layer is not stratified anymore and the 13 C 
pocket which forms at the end of the third dredge-up is polluted with the 14 N from layers 
in the partial mixing zone with an initially larger H/ 12 C ratio. Using the mixing law from a 
stellar evolution sequence of an initially rapidly rotating star the resulting neutron exposure 
is too small to allow any production of s-process elements. 

We have constructed synthetic models with a parameterized range of mixing efficiencies 
during the interpulse phase. In these models we find a spread of neutron exposures in the 
s-process layer. Hence, we tentatively propose that the signature of rotationally induced 
mixing in AGB stars might be identified with the observationally-inferred spread of neutron 
exposure in the interpulse. The spread of mixing efficiencies may be related to a spread of 
angular velocity gradients in stars of different initial rotational velocities, masses, or evolu- 
tionary status. However, mixing efficiencies which would correspond to the observed spread 
of neutron exposures are substantially lower than the mixing law in our stellar evolution 
model with rotation. In this scenario this might be attributed to some effect of angular mo- 
mentum redistribution and loss due to a physical process missing in our calculations, such 
as efficient third dredge-up, mass loss and magnetic fields. 

In the framework of the radiative interpulse s-process production site our simulations 
establish a relation between the intershell carbon abundance and the range of mixing effi- 
ciencies during the interpulse phase of AGB stars. Any mixing in the s-process layer during 
the interpulse phase and after the initial formation of the partial mixing zone can only 
reduce the final neutron exposure. A spectrum of slow mixing efficiencies during the in- 
terpulse phase can deliver a spread in neutron exposures in accordance with the observed 
spread of s-process indices if the neutron exposure in the limiting case without mixing is 
larger than the largest observationally required neutron exposure (for example case ST • 2 
in Busso et al. 2001). If the carbon intershell abundance is larger than that obtained in 
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models without hydrodynamic overshooting at the bottom of the He-flash convection zone, 
then the 13 C abundance in the PM zone and thus the neutron exposure will be larger as well. 
As discussed in §5.3 AGB models including this effect are in very good agreement with the 
observed carbon abundances of H-free and hot central stars of planetary nebulae, in which 
the progenitor AGB intershell abundance is believed to be seen at the surface. 

In this paper we have used the nucleosynthesis and the observed features of the s-process 
to constrain stellar models. We have demonstrated how different mixing processes in differ- 
ent locations and at different times impact the current model of s-process nucleosynthesis in 
low-mass AGB stars. In the future we should improve the stellar physics in order to inves- 
tigate which processes determine the mixing at the core-envelope interface and why mixing 
efficiencies in our rotating stellar model are not in agreement with the requirements of the 
s-process. The theory of angular momentum transport and loss may need to be revisited. 
The question of the third dredge-up in rotating AGB stellar models needs additional con- 
sideration, and it is likely that the treatment of \x gradients and their impact on mixing will 
play an important role. 

F. H. appreciates support from D. A. VandenBerg through his Operating Grant from 
the Natural Science and Engineering Research Council of Canada. F. H. also thanks the 
Institute of Nuclear Theory at the University of Washington for its hospitality and the U.S. 
Department of Energy for partial support during the completion of this work. The detailed 
and thoughtful comments by the referee have helped us to improve the presentation of this 
paper. 
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A. Details of the computation of [hs/ls] and s/s vs r 

In §2 we have presented a s-process calculation (Figure 2) to establish the relation be- 
tween the neutron exposure and the s-process indices [hs/ls] and log(s/so). We have used for 
this calculation initial mass fractions for the light neutron poison typically encountered in the 
partial mixing zone, including ( 4 He, 12 C, 16 O, 22 Ne) ~ (0.74,0.23,0.01,0.02). In our calcula- 
tion the temperature and density are chosen constant at T = 9.8 • 10 7 K and p = 2000g/cm 3 . 
For intermediate- and heavy-mass isotopes we have used an initial solar abundance distribu- 
tion. The charged-particle reactions listed in § 3 and neutron-capture reactions on isotopes 
up to Pb are considered in the calculation. We have taken an initial mass fraction of 13 C of 
0.03 and assumed only a residual amount of initial nitrogen (X( 14 N) = 1.2- 10~ 5 ). The initial 
amount of protons was set to zero. With this choice a final neutron exposure of r ~ 0.7 
mbarn -1 is reached when all 13 C is consumed. 

In a partial mixing situation a significant amount of 14 N may be present. Due to different 
assumptions about mixing in the underlying stellar evolution model, the abundance of 12 C 
and 16 may be different from model to model. The larger intershell 12 C abundance in models 
with hydrodynamic overshooting increases the neutron exposure that can be generated in 
the partial mixing zone (Lugaro et al. 2003b). However, we show here that the relation 
between r and the s-process indices does not depend on the abundance distribution of light 
poison, for example the presence of 14 N. We demonstrate this by some test calculations and 
by recalling some basic properties of the s-process as layed out in text books like Clayton 
(1968). 

The neutron exposure is defined as r = P N n v^ dt, where N n is the neutron density and 
v t is the thermal velocity. The neutron density is related to the molar neutron abundance by 
N n = NApV(n) A(n) where all symbols have their usual meaning. Since neutron captures 
are faster than any of the charged particle reactions the neutron density is given by the ratio 
of neutron production by sources and destruction by sinks. The presence of light neutron 
poison (sinks) will depress the neutron density (and eventually the neutron exposure) and 
will thereby determine the amount of neutrons available to enter the s-process path. The 
s-process distribution (as given by the [hs/ls] index) as well as the total production of the 
s-process elements (as given by the index log(s/so)) is then given by the number of neutrons 
that are captured by the s-process seed (usually 56 Fe) and other heavy elements. This number 
depends solely on the neutron density and the abundances and reaction rates of seed nuclei 
and all s-process species themselves. While neutron poison may limit the neutron density 
(and eventually r) they can not directly influence the relation of the integrated neutron 
density (r) and the s-process distribution and total enhancement. 

More specifically, the s-process distribution depends univocally on the number of neu- 
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trons, n 56 , captured by 56 Fe and its progeny. If the abundance of 56 Fe is much larger than the 
abundance of its progeny n 56 is equal to the amount of 56 Fe destroyed by neutron captures: 

F(n 56 ) = -F( 56 Fe) = J F( 56 Fe) Y(n) r 56 , (n7) dt, (Al) 

where r 56) („ 7 ) is the rate of the 56 Fe(n, 7) reaction and Y are molar abundances. If we assume 
13 C(a,n) to be the neutron source and 56 Fe(n, 7) and 14 N(n, p) to be the major neutron 
absorbers, then with the equilibrium condition dN n /dt = the total neutron abundance is 
given by: 

Y(n) = Y ( 13 C) Y ( 4 He) r 13 , (a , w) 
1 ' Y(^)r 14M + Y^Fe)r 56Xn ^ 

where ri^ n ,p) is the rate of the 14 N(n, p) reaction and ri 3i ( Q)n ) is the rate of the 13 C(a,n) 
reaction. The abundance of light neutron poison like 14 N determines the neutron density 
on which the number of neutron captures by 56 Fe and progeny depends (Eq. Al). But the 
light neutron poison do not affect the relationship between Y{ji^) and Y(n), i.e. between 
the s-process distribution and the neutron exposure. 

In Figure 13 we show different cases of the same computation of Figure 2 performed 
with a range of initial 14 N mass fractions. Changing the initial mass fraction of 14 N leaves 
almost completely unchanged the relationship between the s-process distribution and the 
neutron exposure. The difference is in the final r value, which is much lower for a higher 
14 N abundance. In the same way the relation of r with the s-process indices is independent 
on the abundance of other light neutron sinks, and rather general. 

It follows that the parameter r determines univocally the s-process distribution. Note, 
that the conclusions reached in § 7 do not change if instead of r we quantify the s-process 
with the parameter (L/G), which was introduced in §4, and which is a measure of the 
number of neutrons captured by the heavy sink particle 63 G, i.e. is similar to considering the 
number of neutrons captured by 56 Fe. In the synthetic models presented in § 7 for a range 
of log Dip rs ... — 2. we obtain a range of r and correspondingly a range of (L/G) (see 
Table 3). 
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Table 1. Effect of the 13 C-pocket mass on the s-process efficiency and abundance 

distribution 



IP 


BM/10 


BM/2 


BM (a) 


BM-2 


[ls/Fe] 


0.03 


0.13 


0.23 


0.38 


[hs/Fe] 


0.06 


0.25 


0.40 


0.61 


[hs/ls] 


0.03 


0.12 


0.17 


0.23 



( a )BM: benchmark 13 C pocket corresponding to case with hydrodynamic overshoot with 
/ = 0.128. 



Table 2. Mixing coefficients and Lagrangian mixing efficiency in the partial mixing layer 

of the AGB model including rotation 






dm/dr^ 


log/# 


log £> m 




Am mix 




M /cm 


[D T } = cm 2 /s 


[D m ]= M|/yr 


yr 


M 


0.034 


5.4 • 10" 14 


~ 5.0 


-14.0 


1440( c ) 


3.7- 10- 6 


0.232 


9.2 • 10~ 12 


~ 1.5 


-13.1 


8560^ 


2.7 • 10~ 5 


0.557 


3.6 • 10- 11 


~ 2.8 


-10.6 


2900< d > 


2.7 • 10~ 4 



( a ^at m r = 0.745965 M ; ^max, see Figure 8 for the overall distribution in mass; (c)time 
since thermal pulse; (d)nuclear lifetime of 13 C against a-capture r a ( 13 C). 
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Table 3. Maximum r and (L/G) obtained for different log Dip 



log Dip 


T 


(L/G) 


[D T ] = cm 2 /s 


mbarn^ 1 


number ratio 


-0.0 


0.12 


4 


-f.O 


0.33 


18 


-2.0 


0.64 


55 


no mixing 


0.72 


70 
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Boundary of envelope convection 
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Fig. I. — Schematic representation of the s-process in the interpulse phase of TP-AGB stars 
in a space-time diagram. The ordinate covers the mass range between the H- and He-burning 
shells (intershell). The time marks represent a rough estimate. 
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Fig. 2. — s-process indices log(s/s©) and [hs/ls] in s-processed material as function of the 
neutron exposure from a network calculation including all relevant light and heavy elements 
from H to Pb. The initial 13 C mass fractions is 0.03 and neutrons are released by the 13 C(a,n) 
reaction at a constant temperature of 9.8 • 10 7 K. The neutron exposure corresponding to 
the shaded areas can be excluded as the the dominant contribution in the s-process site (see 
text). As discussed in Appendix A the presence of light neutron poison like 14 N would limit 
the neutron exposure that can be achieved but the relations of r and the shown indices are 
only very weakly effected. 
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Fig. 3. — Comparison of the s-process parameters n cap and L/G as functions of r. The 
quantity n cap is computed with a full s-process network calculation while the ratio L/G 
comes from the calculation in which the s-process is approximated with a sink. The first 
increase of n cap with respect to L/G at very low r is due to neutron captures on isotopes 
from 56 Fe to 62 Ni. From r ~ 0.5 mbarn" 1 the results obtained with the neutron sink show 
larger and larger differences from the results obtained with the full s-process network. The 
neutron sink representation is a valid approximation for r < 0.6 mbarn -1 , i.e. in the range 
of interest for solar metallicity stars. 
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Fig. 4. — Neutron density (solid line, left ordinate) and Maxwellian-averaged neutron capture 
cross section (dashed line, right ordinate) of sink particle 63 G (described in §4) as a function 
of time from s-process calculations as those performed for Figure 2 with a network including 
all heavy elements up to Pb. Three test calculations are presented starting with three 
different amount of 13 C mass fractions (indicated by the labels). The variation of the cross 
section of 63 G reflects its dependence on the abundance distribution of the species represented 
by the sink particle. 
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Fig. 5. — The evolution in time of the neutron exposure r for the three calculations shown in 
Figure 4 is represented by the solid lines (initial 13 C mass fractions indicated by the labels). 
Each calculation is compared to three tests computed with different values of the Maxwellian- 
averaged neutron capture cross section of the neutron sink reaction cx( 63 G) described in §4. 
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Fig. 6. — for caption see page 45 
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Fig. 7. — Evolution of angular velocity profiles through the 25 th interpulse-pulse cycle. 
Labels give interpulse phase = (t — to)/ At, where the previous pulse, i.e. the maximum 
He-burning luminosity, has occured at to; the interpulse period is At = 3 • 10 4 yr. The mass 
range covers the intershell region from the top of the C/O-core (location of the He-shell) to 
the bottom of the convective envelope. 
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Fig. 8. — Evolution of mixing coefficient due to rotationally induced mixing during the 
interpulse phase; labels give interpulse phase (f> — (t — t )/At, where the previous pulse, 
i.e. the maximum He-burning luminosity, has occured at t ; the interpulse period is At = 
3 • f0 4 yr- Two different mixing periods resulting from rotation occur in the partial mixing 
zone during the interpulse period (see text). 
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Fig. 6 - Abundance profiles in the partial mixing zone at three times during the 7 th 
interpulse phase of the model with hydrodynamic overshoot and no rotation. Top panel: 
first post-processing model after the end of the third dredge-up phase; middle panel: 10% 
of 13 C is burned by 13 C(a, n) 16 0; bottom panel: neutron exposure r at end of the s-process 
in the partial mixing zone when the 13 C neutron source is exhausted. 

Fig. 9 - Same as Figure 6 for the simulation at the 25 th interpulse phase including 
rotation and no hydrodynamic overshoot. The middle panel shows the same interpulse 
phase as the middle panel in Figure 6 with respect to the nuclear lifetime of 13 C against 
a-capture. 
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Fig. 9. — for caption see page 45 
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Fig. 10. — Maximum neutron exposure in the stellar layer where the s-process occurs during 
the TP-AGB interpulse period as a function of phase 0. The solid line is the result of detailed 
computation with rotational mixing (shown in Figure 9) and the long- and short-dashed lines 
show test calculations (see text, §7). 
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Fig. 11. — Final neutron exposure profiles in the s-process layer of simulations with a 
synthetic mixing law, combining the initial formation of a partial mixing layer by overshoot 
with a constant interpulse mixing coefficient that mimics the effect of shear mixing (see text 
for details). 
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Fig. 12. — Profile of 15 N mass fraction across the partial mixing zone at a time close to the 
maximum neutron density in the simulation with rotation presented in §6. the long-dashed 
and dotted lines clarify the contributions by two reaction chains discussed in the text. For 
comparison the 15 N profile for the same time is shown for the test case rotl presented in § 7. 
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Fig. 13. — Same as Figure 2, but each case with a different initial abundance of 14 N. The 
labels indicate the logarithm of the 14 N mass fraction. Case (a) has been computed with 
intital abundances X( 4 He) = 0.95, X( 12 C) = 0.03, X( 14 N) = 1.2 • 10~ 5 , and p = 20g/cm 3 . 



